Light emission from Na/Cu(lll) induced by a scanning tunneling microscope 
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Measurements of light emission from a scanning tunneling microscope probing a Na overlayer on 
the (111) surface of Cu are reported along with results of a model calculation that essentially agree 
with the experimental ones. The observed light emission spectra show two characteristic features 
depending on the bias voltage. When the bias voltage is smaller than the energy of the second 
quantum well state formed outside the Na overlayer the light emission is due to a plasmon-mediated 
process, while for larger biases light emission is mainly caused by quantum well transitions between 
the two levels. 

PACS numbers: 73.20.At, 68.37.Ef, 73.20.Mf, 73.21.Fg 



I. INTRODUCTION 

The formation of quantum well states (QWS) near the 
Fermi level is a characteristic feature of the electronic 
structure of alkali overlayers on the (111) surfaces of no- 
ble metalsiSiSi^i^ In these systems electrons may be 
confined to a narrow surface region; they cannot escape 
from the surface because of the vacuum barrier, and at 
the same time they cannot propagate into the noble metal 
because the periodic potential in the bulk creates a local 
band gap near the L point of the Brillouin zone. Al- 
kali overlayers thus offer an interesting opportunity to 
study confined electron systems in metals and they have 
attracted considerable interest. 

Up to now, most investigations of alkali overlayers 
have been carried out using photoemission spectroscopy 
(PES), inverse photoemission spectroscopy (IPES) and 
two-photon photoemission (2PPE)i2iLSi2iiS Thus the en- 
ergies of the quantum well states have been studied as 
a function of substrate and overlayer material as well as 
overlayer thickness (or coverage) . In addition to measur- 
ing the position in energy of the quantum well states, the 
lifetimes of these states have also been addressediii*iSii^ 
The quantum well states overlap in energy and space with 
three-dimensional states, producing interesting possibil- 
ities for quasiparticle decay to proceed simultaneously 
through both two- and three-dimensional channels as well 
as electron-phonon scattering. 

In a recent experiment the system Na/Cu(lll) was 
studied with an alternative experimental technique, light 
emission induced by electron injection from the tip of a 
scanning tunneling microscope (STM).^'* This technique 
has previously been used in studies of quite a few systems, 
most notably noble metal surfaces and semiconductors.^^ 
In those cases photons are generated by a fraction of 
the tunneling electrons that undergo inelastic tunneling 
processes in which an amount of energy, limited to eU 
with U the bias voltage, is transferred from the elec- 
tron to the photon. If the substrate is metallic, the 



rate of spontaneous light emission is increased compared 
with the case of inverse photoemission from an isolated 
surface.i6d7,i8d9,20,2i,22,23,24,25 ^his is due to enhanced 

vacuum fluctuations of the electromagnetic field in the 
tunnel gap between tip and sample as a result of the for- 
mation of localized interface-plasmon modes there. For 
semiconductor samples, on the other hand, electromag- 
netic effects are less important, instead the light emission 
is due to interband transitions in the semiconductor 
The STM tip serves to locally inject or generate minority 
charge carriers which then recombine with majority carri- 
ers giving rise to luminescence. Similar studies have also 
been carried out on semiconductor surfaces semi- 
conductor quantum wells52. and quantum dotS)2i and ad- 
sorbed molecules 

Interestingly enough, the light emission from 
Na/Cu(lll) reported in Ref. Q, appears to involve 
both mechanisms described above. For a wide range of 
coverages, Na on Cu(lll) exhibits several QWS. The 
lowest one, QWSl, has an energy close to the Fermi 
level (see Table QJ. At the coverages to be discussed 
below, it is either occupied or unoccupied. The other 
state, QWS2, lies well above the Fermi energy. As long 
as the bias voltage is low enough that electrons cannot 
be injected into QWS2 light emission is still possible, 
and it proceeds through the same "plasmon-mediated" 
mechanism as in the case of clean noble-metal surfaces. 
In this case electrons undergo an inelastic transition 
from a filled state in the tip to the empty QWSl and 
a photon is emitted. This gives a fairly broad peak in 
the light emission spectrum. Once the bias voltage is 
high enough that electrons can be injected directly into 
the upper quantum well state (QWS2) the emission 
mechanism changes. A large part of the tunnel current 
will now go through QWS2, and light emission will 
mainly be due to transitions between QWS2 and QWSl, 
yielding a fairly sharply peaked spectrum. Of course, 
since the QWS wave functions are to a large extent 
confined to the region of space between the tip and 
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sample, the light emission rate is still enhanced by the 
electromagnetic fluctuations there. 

In this paper we will support the scenario outlined 
above by model calculations that lead to very good qual- 
itative agreement between experiment and theory. We 
will also present additional experimental results. 

The rest of the paper is organized in the following way. 
In Section^ we outline the experimental setup, and the 
experimental results are presented in Sec. IIIII Section 
IIVI gives a comprehensive description of the model cal- 
culation that we have carried out. Next, in Sec. we 
compare theoretical and experimental results and discuss 
the implications. The paper is summarized in Sec. IVII 
Finally, an Appendix presents a derivation of the rate 
of light emission based on the Keldysh non-equilibrium 
Green's function technique. 



II. EXPERIMENTAL DETAILS 

The experiments were performed with an ultra-high 
vacuum (UHV) STM operated at a temperature T = 
4.6 Ki^ Photons in the energy range 1.2 eV <\\v < 3.5 
eV were detected with a lens-system in UHV coupling 
the light to a grating spectrometer and a liquid nitrogen 
cooled charge-coupled device camerai^ All spectra have 
been corrected for the wavelength dependent detection 
efficiency. Throughout the measurements it was verified 
that the surface structure was not modified during data 
acquisition. 

W tips were prepared by electrochemical etching 
and subsequent sputtering and annealing in UHV. The 
Cu(lll) surface was cleaned by repeated cycles of Ar-ion 
bombardment and annealing. Na films were evaporated 
from outgassed SAES Getters sources onto the Cu crystal 
held at room temperature. A quartz crystal microbalance 
was used to estimate the coverage @ which was further 
calibrated by the known binding energies of the quantum 
well states iSiLSiSiili After preparation at room tempera- 
ture the sample was transferred to the STM and cooled 
to T = 4.6 K. 

Following Ref. we define 1 monolayer (ML) as the 
most densely packed structure of the first Na layer, 
namely a (3/2 x 3/2) mesh. This pattern corresponds 
to 4 Na atoms per 9 first layer Cu atoms. 



III. EXPERIMENTAL RESULTS 

Figure ^ shows fluorescence spectra recorded over a 
range of sample voltages U from Cu(lll) covered with 
1 monolayer of Na. The emission from these surfaces 
is comprised of two distinct components. First, there is 
an emission band which is most clearly observed at low 
U . Its maximum shifts linearly to higher energies as U 
is increased. This emission is similar to the plasmon- 
mediated emission observed from noble metal surfaces 
due to inelastic tunneling processes>i2iiSii2i2ii21 



A new spectral structure, the position of which only 
weakly depends on the sample voltage ?7, emerges at 
hv ~ 2.4 eV when the bias is raised to U ^ 2.5 V. From 
its intensity, assuming isotropic emission, a quantum effi- 
ciency of roughly 10^^ photons per electron is estimated 
which is higher than typical values for conventional in- 
verse photoelectron spectroscopy. 

Qualitatively similar observations were made at other 
coverages (Fig.|2l. In addition to a plasmon-related emis- 
sion (circles) at low bias voltages we find emission (dots) 
at higher photon energies which are almost independent 
of the applied sample voltages. The positions of these 
emission features vary with the Na coverage: hv k, 1.7, 
2.5, and 2.1 eV at coverages 0.6, 1, and 2 ML, respec- 
tively. Comparing these photon energies to the data of 
Table n the emission is assigned to interband transitions 
between quantum well states. Small deviations between 
the expected transition energies {E2 — Ei) and the mea- 
sured photon energies are in part due to the electric field 
of the tip which causes a Stark shift i^2i^2iia This shift is 
strongest for higher energy states. 

Closer inspection reveals a distinct difference of the 
data from 0.6 and 2 ML compared to the 1 ML case. At 
1 ML coverage, a sharp drop on the high energy side of 
the emission occurs owing to energy conservation given 
by the condition hv — eU. The emission from the 0.6 
and 2 ML Na films exhibits a similarly sharp drop at a 
different energy, hv ~ eU — Ei, where Ei is the energy 
of quantum well state 1 relative to the Fermi level. 



TABLE I: Energies of quantum well states of Na on Cu(lll) 
in eV relative to the Fermi energy Ef as obtained from tun- 
neling spectroscopy of the differential conductance. At a cov- 
erage of 1 monolayer, the lowest state is occupied {Ei < 0). 
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IV. MODEL AND THEORY 
A. General considerations 

In this section we will use a rather simple, basically 
one-dimensional model of the W-tip/Na/Cu(lll) system 
that nevertheless captures the essential physics of the 
experiment and yields calculated spectra that can explain 
the experimental results. 

The general framework of this calculation is based 
upon what we have used in earlier calculationsiS*^ of 
light emission from noble metal surfaces. As the more 
detailed derivation in the Appendix shows, the intensity 
of the emitted light (per unit photon energy and solid 
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FIG. 1: Light emission spectra of 1 ML Na measured for a 
series of sample voltages between L815 V and 2.815 V with 
a tunneling current of 10 nA. Small apparent undulations at 
low (< L6 eV) and high (> 2.6 eV) photon energies are due 
to counting statistics. There is a qualitative change once the 
voltage reaches « 2.5 V: Below, the peak emission shifts with 
the sample voltage; above, the emission maximum remains at 
hv « 2.4 eV. 



angle) can be calculated from the expression 
dP 27r ^ , . ,T|G(tj)|2 ^^2 
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Thus, the intensity is found from a summation over filled 
initial electron states i and empty final states /, and jji 
is the current matrix element between these states, 
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In Eq. eo is the permittivity of vacuum and c the 
speed of light, and G, finally, is an electromagnetic en- 
hancement factor. Note that Eq. agrees with Eqs. (6) 
and (7) of Ref. '4l', considering that the formulas in that 
work employed CGS units. 

Since the calculation of G(w) in Eq. Q has been 
treated in earlier papers (see, for example, Ref. l42f) we 
will just give a brief outline of the procedure here. As 
discussed in the Appendix, a reciprocity relation makes 
it possible to interchange the source and detection points 
in the electromagnetic calculation. Thus, since we are in- 
terested in evaluating the light emission intensity found 
at an observation angle 9—\ rad., we let an incident elec- 
tromagnetic wave hit the tip-sample system from exactly 
that direction and calculate the electric field enhance- 
ment at a point between the tip apex and the surface, 
where it reaches its highest values. The field in the cavity 
between the tip and sample is essentially constant along 
the direction normal to the surface, however, in the lat- 
eral direction the field enhancement begins to drop off 
at distances exceeding 2-3 nm from the symmetry axis. 
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FIG. 2: Characteristic spectra, obtained with a low sample 
voltage (circles) for which injection of electrons into the upper 
quantum well state is not possible, and a high sample voltage 
(dots) where electrons can be injected into the upper QWS. 
The three panels display results from samples covered by 0.6, 
1, and 2 monolayers of Na. A small dip of the intensity at 
~ 1.3 eV is due to a sharp absorption of the optical fiber used 
in the experiment, which is not fully corrected for. 



The optical properties of tip and sample are modeled us- 
ing macroscopic dielectric functions for Cu and W found 
from Ref. (ij- We have approximated the optical prop- 
erties of the Na layer by the Cu dielectric function since 
the treatment of such a thin layer in terms of macroscopic 
dielectric functions would not be very reliable44 The tip 
is represented by a sphere with radius i? = 30 nm. 



B. Model potential 

We use a one-dimensional model potential, illustrated 
in Fig. 131 to calculate both the tunneling current and the 
matrix elements that set the light emission rate. The 
parameter values that enter this model originate from 
the physical properties of bulk Cu and bulk Na that are 
most relevant to the problem at hand. Similar models 
have been used in earlier calculations^^i^ (see also Ref. 

S3). 

Of course, using a strictly one-dimensional model cor- 
responds to a situation where both electrodes, sample 
and tip, are completely flat. This is not the real situa- 
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tion in an STM experiment. Thus, we use the ID model 
to calculate a current density which then multiplied by a 
suitable effective area yields the tunnel current. A similar 
procedure is applied for the matrix-element calculations 
as well. 

Inside the copper sample the potential is modulated 
along the (111) direction (normal to the surface) as 

Vcu{z) = Vin(e's^ + e-^s^) = 2 cos (52), (3) 

where the reciprocal lattice wave vector g = 27r/acu, and 
acu = 2.08 A is the inter-plane distance in Cu in the 
(111) direction. The top layer is centered at z = 0, and 
the copper sample is assumed to end at z = zq^ = acu/2. 
The value of Vm is taken to be 2.5 eV.4£ The potential 
in yields a band gap for energies (related to electron 
motion in the z direction) between Eg — Vm and Eg + 
Viii, where Eg = n^{g/2f/{2m) « 8.7 eV is the kinetic 
energy of a free electron with a wave vector at the L 
point on the Brillouin zone boundary. Thus, just as in 
copper our model has a 5 eV wide band gap at the L 
point. This feature is crucial in forming quantum well 
states in the overlayer system. Experimentally the energy 
difference between the lower band edge and the Fermi 
level is 0.9 eV, we therefore put the Cu Fermi level at 
EF = Eg- Vm + 0.9 eV = 7.1 eV. 

We treat the sodium overlayer as a bulk, free-electron 
metal. The thickness is assumed to correspond to 2 
monolayers (6.13 A)4& The constant value of the poten- 
tial in the sodium is set to Vrja = 3.9 eV. The choice 
of this value is based on the fact that the Fermi energy 
of bulk Na is 3.2 eV, so that now an electron with the 
kinetic energy of 3.2 eV in the Na layer will have a total 
energy of (3.9 + 3.2) eV = 7.1 eV, i.e. identical to Ep- 

Moreover, in order to get a tunnel current from elec- 
trons with energies in the Cu band gap, it is necessary to 
add a negative imaginary part to the potential. We have 
chosen to put this in the sodium layer, thus the full Na 
potential is 

^Na-ir, with r = 0.1eV. 

r in an approximate way represents the effects of 
electron-phonon, electron-electron, and interface scatter- 
ing that in the real system eventually remove electrons 
from the quantum-well states near the surface. If the po- 
tential had this value in all of space, the electron density 
due to a particular state would decay in time as e"^^*/'^. 
Thus the lifetime would be r = ?i/(2r) « 3 fs, and the 
peak in the spectral function associated with the state 
would have a full width at half maximum of 2r. It must 
be kept in mind that an electron in a QWS only spends 
part of the time (sa 50 % for QWSl) in the Na layer. Our 
choice for F would thus give a lifetime of some 6 fs, and 
a hue width of 0.1 eV for QWSl. These values are com- 
parable with measured values found in the literaturei^ 
We have for simplicity assumed that the potential is 
constant in the W tip. This is of course a rather crude 
approximation, but not a crucial one since the electronic 



structure of the tip is not of primary importance for 
the analysis of the experiment at hand. With the free- 
electron band width W = 8.0 eV, the potential in the 
tungsten tip at zero bias is Vtip = Ep — W = —0.9 eV, 
however, when the tip is biased the potential is given by 
Vtip = EF~W + eU. 

Finally, we need a barrier potential to use in vacuum 
between the Na overlayer and the tip. This potential es- 
sentially consists of two parts: a tilted square barrier and 
image potential contributions. The tilted square barrier 
can be written 

^tiit(^) = '-^^^{EF + 4>n,) + ^—^{EF + eU + 4>tip), 

(4) 

where ztip is the coordinate of the tip apex, ZNa is the 
coordinate of the Na overlayer surface, and d is the tip- 
sample separation. Here 0Na = 2.7 eV and 0tip = 5.2 eV 
denote the Na and tip work functions, respectively, and 
U is the bias voltage. The image contributions are 

" "47reo 4(z-z;j™) " 47reo 4 (z™ - z) " 
We have chosen the image plane position as 

4" - ^Na = -eV[16^eo('^Na + Ef~ VNa)] (6) 

in the sample. This guarantees that the total potential 
Vtiit + Vim equals Vwaj the potential inside Na, at z = z^a 
provided the tip is far away. In the tip, following the same 
reasoning, 

zj™ - ztip = eV[167reo((/.tip + W)]. (7) 

The resulting potential is illustrated in Fig. |21 for two 
different values of the bias voltage. 

The model potential described above yields an elec- 
tronic structure of the Na overlayer that is in reasonable 
agreement with experimental observations. In particular, 
using this model we are able to reproduce the essential 
features found in the experiment by Hoffmann, Kliewer, 
and Berndt.i" At the same time it should be pointed out 
that the model is not detailed enough to reproduce ex- 
actly the same energy level positions as found experimen- 
tally. Given that the Na layer is modeled using bulk pa- 
rameters, the theoretical results should be more accurate 
for thicker overlayers. We have therefore concentrated on 
calculating results for the 2-monolayer case. 

C. Wave functions 

To calculate the tunnel current and subsequently the 
light emission intensity we must solve the Schrodinger 
equation for the electron wave functions. Since we have 
chosen to work with a potential that is translationally 
invariant in the directions perpendicular to the tunneling 
direction, we can write all wave functions on the form 

V'(r) = V(^)e'^"-*'". (8) 
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FIG. 3: Illustration of the model potential at two different 
values of the bias voltage. In addition the wave functions cor- 
responding to the two quantum well states have been calcu- 
lated using boundary conditions corresponding to an electron 
impinging on the barrier from the tip side when [/ = 2 V. 
Note that the wave function at the QWS2 energy penetrates 
the barrier with high probability, however, since this energy 
here lies above the tip Fermi energy, this does not contribute 
to the tunnel current. 



In the copper sample, in view of the potential given by 
Eq. we make the nearly- free electron model Ansatz 
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(3e 



i{k-g)z 



(9) 



for the wave function. Inserting this into the Schrodinger 
equation one finds that there is a band gap for energies 
(related to the motion in the z direction) between Eg — 
Viii and Eg -\- Vm. In this energy interval there are 
no bulk states in the copper, however, one can still have 
surface states for which k = p ~ iq with p = g/2, and 



2m{JAeEg + V^\^ -Eg- s)/h. 



Thus, the wave function envelope decays exponentially 



(with fctip = y/2m{e — Vtip)/h) at z = ztip. The coeffi- 
cients C and F, as well as A and B in Eq. ifTIUl . must be 
determined by wave function matching. 



D. Tunnel current 

Before we can evaluate the matrix elements entering 
Eq. ^ we have to calculate the tunnel current Jdc, or 
rather, determine the tip-sample separation d that yields 
a certain, set value for /(jc- The probability current den- 
sity in the z direction associated with one particular wave 
function can be written 

..^^^Re[^*(.)|^(.)]. (12) 

Here the first two factors serve as normalization. A^s is 
the effective area of the tunnel contact, which we also 
use as a normalization area for the wave functions!^ L 
is a normalization length (in the tip), and by dividing 
by \F\'^ wc normalize the current to the current carried 
by the wave impinging on the barrier from the tip side. 
For a real-valued potential the probability current is 
independent of the z coordinate. In the present case, 
with a potential that has a non-zero imaginary part in 
the sodium layer, the probability current is independent 
of z to the right (in the barrier and tip) of the Na layer, 
and can be evaluated anywhere in that part of space. 

vanishes, on the other hand, in the copper since we 
cannot have any propagating states in the energy gap 
there. A tunnel current is flowing across the tunnel gap 
only because of the scattering processes that eventually 
scatter electrons out of the quantum well states confined 
to the Na overlayer and surface region of the copper. 

The total electric current is obtained by summing over 
all contributing states to get a total current density and 
then multiply by the electron charge and the effective 
tunneling area 



idc = -e^cff 



Jtp- 



(13) 



qz 



and the electron is confined to the region near This sum can then be turned into an energy integral. 



as e 

the copper surface. Since we use a positive value for the 
corrugation of the potential (Vm = 2.5 eV) states near 
the bottom of the gap will be p-like, whereas states near 
the top of the gap are s-likei^ 

The wave function in Na takes the form 



ipiz) = Ae'''''-^ + Be 



(10) 



where ^Na = \/2m(e — V^a. + iT)/h and the branch cut 
of the square root function is placed along the negative 
real axis so that the imaginary part of /cNa is positive. In 
the barrier region, the wave function must be integrated 
numerically and then joined to the tip wave function 



V'(z) = Ce^'^'^p^ + Fe 



-iktipZ 



(11) 



E. Intensity of emitted light 

To calculate the differential power of the emitted light 
we must carry out the sum over initial and final electron 
states as indicated in Eq. In contrast to the case of 
light emission from noble metals the final electron state 
is in this case, at least in one sense, discrete. Light is pri- 
marily emitted while the electron traverses the vacuum 
barrier and then QWSl is in practice the only possible 
final state (assuming parallel momentum is conserved); 
it dominates the local density of states completely in the 
energy range just above the substrate Fermi level. 
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Therefore in the calculations, we have solved for a 
bound state to represent ipf. This state is calculated 
with essentially the same potential as was used for the 
tunnel current calculation, however, in this case we set 
r equal to zero, and determine an energy eigenvalue by 
requiring that the wave function is decaying well inside 
the tunnel barrier. In the x and y directions this bound 
state wave function is assumed to behave like a plane- 
wave state with a certain momentum. Furthermore, in 
the following we will assume that this state is almost al- 
ways unoccupied so that tunneling electrons can make 



inelastic transitions into it and emit light at the same 
time. This last assumption is reasonable because even if 
most of the tunnel current passes via QWSl, the lifetime 
of this state is of the order 10~^^ s while with a tunnel 
current of 10 nA the delay between each tunnel event is 
1.6 X 10"" s. 

In our one-dimensional model the momentum parallel 
to the interfaces and the electron spin is conserved in the 
inelastic tunneling process. The sum over initial and final 
states in Eq. ^ then reduces to integrals over parallel 
and perpendicular momentum. 



dP 



d{fujj)(Kl 



uj''\G{uj)\ 
87r2eoc3 
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"27 



\jM'6iE,~Ef~hu). (14) 



Here the integration over parallel momentum can be turned into an energy integration over an interval starting at zero 
and ending at the maximum energy E'n.max that the electron can have due to the motion parallel to the interfaces. 
This energy is given by the tip Fermi energy, i.e. Ep + eU minus the z motion energy Ef + huj in the initial state, i.e. 



El 



I, max = Ep + eU — Ef — huj. We then get 

dP ,C<j2|G(w)l2 



= 2- 

d{huj)dQ Sn^eoc 



mE\ 



II ,max 



0(^^11, 



(15) 



where Q denotes a step function. It remains to carry out the kz integration, and also this is, in view of the 5 function in 
the integrand, straightforward. If we let Ef, stand for the band bottom energy in the tip, one has Ti^ hi /{2m) ^E^-Eb. 
It is then possible to show that 



S{E^ - Ef - Tioj) 



5\kz + \ 2m{Ef + nuo- Et)/h' 



^^2m{Ef + riLo-Eb) 
and when this is inserted into Eq. H15|l we get the final result 



dP 



d{huj)dfl 



A, 



e^uj^\G{Lu)\^ {Ef + eU ~ Ei - hLj)Q{EF + eU - Ei - tnu) 



647r4en?ic3 



^2to(£'i + nuj- Eb) 



dz 



dz ■' dz 



(16) 



(17) 



Note that we have made the replacement Ef Ei in 
the last equation because the final state here is identi- 
cal to the lowest QWS. The z integration is limited to 
the vacuum part of space, since the electromagnetic field 
enhancement is much higher there than in the tip and 
sample. 



V. RESULTS AND DISCUSSION 

Figure 01 shows spectra calculated from the model of 
Sec. IIVI for a number of different bias voltages. The cal- 
culated tunnel current was kept constant at Idc — 10 nA 
in all cases by changing the tip-sample separation. 

We see that these spectra have a number of general 
qualitative features in common with the experimental 



spectra shown in Figs.n]and[21 At first, for bias voltages 
< 2.5 V, the spectrum has a relatively broad peak that is 
cut off on the high-energy side. Since QWSl is the final 
state in the light emission process the maximum photon 
energy equals eU — {Ei — Ep), and not eU, [cf. Eq. ((T7|) ]. 

Once the bias voltage becomes high enough that elec- 
trons can be injected into the upper QWS most of the 
tunnel current will take this path. At the same time the 
tip-sample separation increases to maintain the tunnel 
current at a constant value. Quite naturally the electrons 
injected into the upper QWS will also dominate the light- 
emission processes at these higher voltages. Most of the 
emitted photons will have an energy closely correspond- 
ing to the energy difference between the two quantum 
well states. 

Returning to the model, these facts can be understood 
as follows. As long as the tip Fermi energy lies below 
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FIG. 4: Calculated light emission spectra for the Na/Cu(lll) 
system with an overlayer thickness corresponding to 2 mono- 
layers. The spectra were calculated for a series of bias volt- 
ages, while the tunnel current was kept fixed at 10 nA by 
varying the distance between tip and sample. 



QWS2 all initial state wave functions will decrease expo- 
nentially upon traversing the tunnel barrier. Therefore 
wave functions "02 , corresponding to a broad range of en- 
ergies give comparable contributions to the integral in 
Eq. 117|) . The shape of the spectrum in this case is mainly 
determined by the phase space factor {Ep+eU — Ei — huj) 
and the field enhancement factor |G(w)|^. 

For higher bias voltages, when Ep + ell exceeds E2, it 
is instead the last factor in Eq. H17I) , the matrix element 
integral, that determines the spectral shape. The final 
state remains the same and therefore, to have a large 
matrix element, the initial-state wave function ipi must 
be large in the part of space where the final state resides. 
This will happen when the initial state energy coincides 
with E2 since then resonant tunneling into QWS2 is pos- 
sible and tpi will look like the upper wave function illus- 
trated in Fig. 121 Even if there are plenty of initial states 
with energies both somewhat below and above E2, they 
will not at all give as large contributions to the emit- 
ted light intensity because their wave functions are much 
smaller near the Na overlayer surface. As a result of 
this the light emission spectrum develops a peak around 
hv = huj = E2 — El. 

The overall intensity of the emitted light is smaller at 
the higher voltages. The main reason for this is that 
raising the voltage moves the peak of the light emis- 
sion spectrum away from the frequency range near 2 
eV where a tungsten tip and a copper sample have an 
interface-plasmon resonance causing resonantly enhanced 
light emission. At a photon energy of 2.3-2.4 eV there 
is still a considerable field enhancement between the tip 
and sample, but G is typically down by a factor of 2 com- 
pared with the resonant case. Thus, the peak in the light 



emission spectrum at hv w 2.3 eV is entirely due to the 
special electron structure of the Na/Cu(lll) surface. 

The experimental results also illustrate the last point 
in an interesting way. With 1 ML coverage [Figs.n^^iid 
|2Ib)] the emission peak due to quantum well (or inter- 
band) transitions falls at a high photon energy where the 
electromagnetic enhancement is relatively small. Conse- 
quently, comparing spectra taken at different voltages, 
the plasmon-mediated light emission yields the more in- 
tense peaks in these spectra. However, for 0.6 ML cov- 
erage [Fig. |2Ia)] the quantum well transition occurs at a 
lower photon energy, near the maximum of the electro- 
magnetic enhancement. In this case the "quantum- well 
peak" is more intense than the peak occurring at the 
lower bias voltage. 

In this context let us discuss why it is reasonable to cal- 
culate the electromagnetic enhancement at a point just 
below the tip apex both for plasmon-mediated and QWS 
transition light emission. In the first case, this is natural 
since the light emission event must take place while the 
electron tunnels from tip to sample. In the latter case, 
the electron is injected into a fairly long-lived state, and 
in principle it may end up in a point quite far away from 
the tip apex before a photon is emitted. However, with 
a lifetime of 10 fs and a lateral electron velocity of 10^ 
m/s the distance traveled by the electron is no more than 
1 nm, thus it would still be in a region where the elec- 
tromagnetic enhancement has not dropped off very much 
(cf. the discussion in subsection IIV All . 

From the calculated results, one can estimate the quan- 
tum efficiency of the light emission process to be 10~^ 
emitted photons per tunneling electron. This number 
compares well with the experimental result. It is consid- 
erably larger than what is observed for inverse photoe- 
mission processes at a single surface, however, the quan- 
tum efficiency of STM-induced light emission from most 
notably Ag samples may reach values between 10^"* and 
10-3. 



VI. SUMMARY 

In summary we have interpreted experimental obser- 
vations of STM-induced light emission from the quan- 
tum well system Na on Cu(lll) using model calculations 
of the electronic structure and the optical properties of 
the tip-sample region. The main features of the exper- 
imental data, namely two distinct spectral structures, 
and their intensity variation with the tip-sample volt- 
age are reproduced by the calculations. The emission 
lines are attributed to (a) emission from a localized plas- 
mon which is excited by inelastic tunneling to a quantum 
well state and (b) transitions between two quantum well 
states. The electromagnetic enhancement present in the 
tip-sample cavity substantially enhances the intensity of 
the emission and explains the observed, high quantum 
efficiencies. 
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APPENDIX 

In this appendix we derive Eq. |^ using the Keldysh 
non-equilibrium Green's function (GF) technique 
We wish to calculate the intensity of the spontaneously 
emitted light at a detection point rg far away from the 
STM tip and sample. Since all photons appearing at this 
point that are of interest to us originate from the region 
around the STM tip, the intensity can be found by multi- 
plying the electromagnetic energy density, excluding the 
zero-point energy, by the speed of light. The energy den- 
sity can be expressed in terms of Keldysh GF's as 

w = «eo / — [L'|£(ro,ro,+cj) + i:)^£;(ro,ro,-cj)], 
Jo 27r 

(A.l) 

where £'^^(ro, Tq, w) and £)^^(ro, Tq, are the Fourier 
transforms of 

D<E{ro,rQ,t) = ~i{Ee{ro,0)Eg{ro,t)) (A.2) 



and 



(A.3) 



respectively. In writing these expressions we have as- 
sumed that the tunnel current causing light emission 
flows in the z direction near the origin and sends out 
p-polarized light with an electric field E = OEg, point- 
ing in the 9 direction in the far field. Consequently, the 
radiated differential power can be written 



dP 



(A.4) 

The two Keldysh GF's in Eq. (|A.4|I actually yield iden- 
tical contributions and in the following we will only deal 
with D<. 

We need to find an expression for the contributions to 
that result from interactions between the electro- 
magnetic field and the electron system in the STM tip 
and sample. This interaction is to lowest order 



H' = TT- y'[A(r„) • p„ + p„ • A(r„)], 
Zm ^ — ' 



(A.5) 



where A is the electromagnetic vector potential (E = 
—dA/dt), the sum runs over the electrons, and r„ and p„ 



are the electron coordinates and momentum operators, 
respectively. By performing an S-matrix expansion of 
to second order in H' we find, using the rules for 
"analytic continuation",^^ that D"^^ can be expressed as 
a product of two photon Green's functions and a current- 
current Green's function, 

x^ba(i"o, i-i, cj)n< (ri, r2, w)D;^£;(r2, Tq, cj). 

(A.6) 

These Green's functions are the Fourier transforms of a 
retarded photon GF 



i?^^(r,r',t) = ^te{t){[Eeiv,t),A,ir',0)]) 
an advanced photon GF 

D-Xsir, r', t) = i9{^t){[A,{r, t),Ee{v\ 0)]) 
and the current-current GF 

n<(r,r',i) = -iO.(r',0)j,(r,t)). 
The current density operator 
-ieTi 



jzir) 



2m 



dz 



dz 



(A.7) 



(A.8) 



(A.9) 



(A.IO) 



where ^ and ^'^ are electron annihilation and creation 
operators.-^ 

At this point we can make a number of approxima- 
tions and simplifications. The inelastic tunneling events 
occur in a very small part of space in the tunnel gap in 
a region of (sub-)nanometer size. The photon Green's 
functions do not vary very much on this length scale, so 
and D"^ can be taken outside the integral, and ri 
and Y2 can be replaced by a fixed point in the tunnel 
gap in these functions. Fourier transformation turns time 
derivatives into frequency multiplications. Therefore the 
photon Green's functions only need to involve the vec- 
tor potential. Moreover, in Fourier space retarded and 
advanced GF's are each others complex conjugates, and 
finally the reciprocity theorem of electrodynamics allows 
us to interchange the source and field points in the re- 
tarded photon GF. This gives us 

D'ea (i"o , i-s , ^^)Dae (r^ , ro , w) = I Dig (r^ , Tq , p . 

(A.ll) 

In a case where parts of space is filled with materials 
characterized by a relative dielectric function er(r) the 
(tensor) photon Green's function D^p solvea^Si^ 



V X V X - er(r)- 



and yields the a component of the vector potential in the 
point r if there is a (5 function current source pointing in 
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the [3 direction in r'. Thus, -D^^(r, r',w) can be calcu- 
lated within the framework of classical electrodynamics. 
The z6 element of interest here can be written 



Dl 



h e 



(A.13) 



where G is an enhancement amplitude. G is given by 
the z component of the local electric field in the tunnel 
gap when a plane wave of unit amplitude incident from 
the direction of Tq impinges on the tip-sample system. In 
free space one would simply have G{9,uj) = s'md. The 
detailed scheme for calculating G has been described ear- 
lier, for example in Ref. 

It remains to evaluate the current-current Green's 
function n^(ri,r2,w) and to carry out the integrations 
over the coordinates ri and r2 in Eq. ljA.61) . A straight- 
forward evaluation yields the result 



J d\i j dV2n<(ri,r2,o^) 



= -27ri?i^ + Ef ~ Ei) (A.14) 

where the sum runs over filled initial electron states i 
with energy Ei and wave function ifii and empty final 
electron states / with energy Ef and wave function ipf, 
and 



—ieh 
2m 



dz 



(A.15) 



By inserting the results of Eqs. (|A.6|) . (jA.llI) . (|A.13|1 . 
and ljA.14|) in Eq. (|A.4|I we arrive at the final result 



dP _ LU^\G{9o,u)\^ 
d{huj)dQ ^ 87r2eoc3 

which is identical to Eq. 



Y,\jM^S{hco + Ef^E,), 

fi 

(A.16) 
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